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ABSTRACT 

Using new Keck DEIMOS spectroscopy, we examine the origin of the steep number counts of ultra- 
faint emission-hne galaxies recently reported by Dressier et al. (2011). We confirm six Lya emitters 
(LAEs) , three of which have significant asymmetric line profiles with prominent wings extending 300- 
400 km s~^ redward of the peak emission. With these six LAEs, we revise our previous estimate of the 
number of faint LAEs in the Dressier et al. survey. Combining these data with the density of bright 
LAEs in the Cosmic Origins Survey and Subaru Deep Field provides the best constraints to date 
on the redshift 5.7 LAE luminosity function (LF). Schechter function parameters, (j)* = 4.5 x 10"'' 
M-pc-^, L* = 9.1 X lO''^ erg s-\ and a = -1.70, are estimated using a maximum likelihood technique 
with a model for slit losses. To place this result in the context of the UV-selected galaxy population, 
we investigate how various parameterizations of the Lya equivalent width distribution, along with 
the measured UV-continuum LF, affect shape and normalization of the Lya LF. The nominal model, 
which uses z ~ 6 equivalent widths from the literature, falls short of the observed space density of 
LAEs at the bright end, possibly indicating a need for higher equivalent widths. This parameterization 
of the equivalent width distribution implies that as many as 50% of our faintest LAEs should have 
Muv > —18.0, rendering them undetectable in even the deepest Hubble Space Telescope surveys at 
this redshift. Hence, ultra-deep emission-line surveys find some of the faintest galaxies ever observed 
at the end of the rcionization epoch. Such faint galaxies likely enrich the intergalactic medium with 
metals and maintain its ionized state. Observations of these objects provide a glimpse of the building 
blocks of present-day galaxies at an early time. 

Subject headings: galaxies: high-redshift - galaxies: evolution - galaxies: formation 



1. INTRODUCTION 

Observations of galaxies in the post-reionization epoch 
(z ~ 4 — 6) measure the growth of young objects in an 
early Universe. Over these redshifts, galaxies assemble 
into the active systems that we observe during the peak 
of cosmic star-formation at z ^ 2 — 3. These galaxies 
drive winds, which play a critical role in the metal enrich- 
ment of t he intergalact ic medium (IGM; [Oppenheimer et 
aL]|2009[ [Martin et al.||2010[ ). Additionally, observations 



of galaxies in the post-reionization era provide a baseline 
for understanding the earliest star formation and rcion- 
ization of the IGM at z > 7. 

Studies of both Lyman Break Galaxies (LBGs; 
Bouwens et al."2007') and Lya emitt ers (LAEs; Rhoads 
et al. 2000; Shimasaku et al 2006; Gronwall et al.||i^007 ' 
Huchi et a'"^'" ' ' 



2010; H u et al.||2010[ [Kashikawa et 

aT 201ip in the post-reionization era provide comple 



mentary views of early star formation and galaxy evo- 
lution. The escape of Lya photons from the interstellar 

^ Some of the data presented herein were obtained at the 
W.M. Keck Observatory, which is operated as a scientific 
partnership among the Cahfornia Institute of Technology, the 
University of California and the National Aeronautics and Space 
Administration. The Observatory was made possible by the 
generous financial support of the W.M. Keck Foundation. 

^ This paper includes data gathered with the 6.5 meter 
Magellan Telescopes located at Las Campanas Observatory, 
Chile. 

^ Department of Physics, University of California, Santa 
Barbara, CA 93106; ahenry@physics.ucsb.edu 

^ Carnegie Observatories, 813 Santa Barbara Street, 
Pasadena, CA 91101 

^ Department of Astronomy and Physics, Saint Mary's 
University, Halifax, NS B3H 3C3, Canada 



media of galaxies is strongly influenced by galactic out- 
flows. Typically, observable emission results when pho- 
tons are redshifted out of re sonance with the systemic 



velocity of the host galaxy (Shapley et al. 2003 Stei 



del et al.||2010[ l2011|. As a result, the comparison be 
tween lAEs and LBGs contains information about star 
formation driven galactic winds. For example, the space 
density of UV-luminous galaxies declines with increasing 
redshift at z ~ 4 — 6, while Lya- luminous gala xies main- 
tain an approximately constan t space density (Bouwens 



et al. 2007 Ouchi et al 



20081. The direct implication 



of these results is that Lya photons more readily es 
cape from galaxies at higher redshifts, possibly due to 
decreasing dust or the increased presence of gaseous out- 
flows. Nevertheless, current conclusions must be drawn 
from limited data, as narrowband imaging observations 
leave the faint-end of the Lya luminosity function (LF) 
unconstrained. 

Looking to the end of rcionization, it is still debated 
whether galaxies produced enough ionizing photons to 
prevent recombination in the IGM. It is often argued that 
galaxies fainter than those present in the deepest Hub- 
ble Space Telescope (HST) surveys are needed to produce 
enough ionizing photons to ba lance recombinations in the 



IGM (|Yan fc Wi^^^^^ 2004 ; 'Sal vaterra et al.|2011 



lator et a l.|20 1ip. Because of this requirement tor fainter 



Fin- 



objects, LAE surveys at this redshift are useful, as they 
can identify galaxies that cannot be detected by other 
means. As an example, an LAE with L = 10'*^ erg s~^ 
and rest-frame equivalent width Wq > 80 A would imply 
a UV-luminosity of Muv > —18.0 — a level too faint to 



be detected in the Hubble Ultra Deep Field ( Bouwens et 



2 



al.|2006 )■ Identifying very faint LAEs can shed new hght 
on the ionizing photon budget by including sources with 
continuum luminosities that have never been probed in 
the epoch of reionization. 

In addition, studies of LAEs in the post-reionization 
era are crucial if we aim to determine the neutral hy- 
drogen fraction from Lya emission. Measurements of 
Lya equivalent widths, line profiles, clustering, and the 
LF will all be impacted by neutral hydroge n (Miralda- 
Escude et al."2000 1 [Furlanetto et aL[|2006| [Pijkstra et] 
al. 2007a b; McQumiret al.| |2007| |Dayal et al.| |2008[ 
Sclienker et aT 201ip , but quantifying these effects will 
require robust measurements in the post-reionization era 
for comparison. Because galaxies in a partially neutral 
IGM are thought to reside in bubbles ionized by their 
own sta r-formation, lumin osity dependent trends are ex- 
pected (Ono et al. 2011). By measuring only bright 
LAEs, narrowband imaging provides an incomplete view 
of these star-forming galaxies. 

Clearly, samples of faint LAEs are required to better 
understand galaxies and the IGM at z > 4. To detect 
these objects, spectroscopic searches are necessary to in- 
crease the contrast with the sky background. Blind spec- 
troscopic searches for LAEs have now successfully iden- 
tified galax ies at 3 < z < 6 ( Santos et al.|2004 Martin 

Sawick 



et al. 
et al. 



2008 



2008 



hereafter Paper I, |Kauch et al. 



Cassata et al.|201ip . In Dressier et al. (2011 



hereafter Paper llj we presented first results from a new 
Multislit Narrowband Spectroscopic (MNS) survey for 
LAEs aX z — 5.7. In this search, we reached unprece- 
dented depths and detected 215 LAE candidates. The 
faintest objects, if they are confirmed, will have lumi- 
nosities that have o nly been seen (at t his redshift) in a 
few lensed galaxies |Santos et aL]|2004[ ). Only 31 (14%) 
the LAE candidates m Paper 11 have Muxes bright enough 
to be selected through narrowband imaging (1.0 x 10^"'^^ 
erg s~^ cm~^^ 

In Paper II, we presented evidence that a high fraction 
of the LAE candidates are indeed LAEs, but spectro- 
scopic confirmation is still required. In the present pa- 
per, we analyze the first results from our spectroscopic 
followup campaign. With a new sample of confirmed 
LAEs in hand, we provide the first measurement of the 
faint-end slope of the Lya LF at redshift 5.7. We demon- 
strate that slit-losses can be robustly accounted for in 
blind-spectroscopic data by using Monte Carlo simula- 
tions within the framework of t he m aximum likelihood 
method. In fact, as we show in 36.21 this method is di- 



rectly applicable to narrowband imaging data, where un- 
certainties often result from a non-uniform narrowband 
throughput. 

This paper is organized as follows: in S|2] we give a 
brief overview of the survey design. In §3 we present our 
spectroscopic followup observations, and in §4 we discuss 
the newly identified LAEs. These LAEs are used, in §5, 
to infer the number of LAEs that should remain in our 
sample, and in §6 we derive the Lya LF by combining 
our data with those fro m the Subar u Deep Field (SDF) 
and COSMOS (Scoville et al.|[2007| ). Prospects for im- 
proving the parameter constraints are discussed in §7. 

® Ground-based narrowb and imaging searches a re limited by 
residual flat-fielding errors ( [Stubbs fc Tonry||2006[ l, so increased 
exposure times will not reveal tainter galaxies. 



The implications of this work, including comparisons to 
other measurements, and a discussion of LAE properties 
can be found in §8. Survey completeness is described in 
Appendix A. Finally, in Appendix B we examine the fore- 
ground (z ^ 0.25— 1.2) emission line galaxies and update 
our model for the interloper fraction that we previously 
presented in Paper II. 

Throughout this paper we use AB magnitudes, and a 
ACDM cosmology: Hq = 70 kms"^ Mpc"\ f^A = 0.7, 
= 0.3. All volumes are given in comoving units, 
unless otherwise stated. 

2. SURVEY DESIGN & FOLLOWUP STRATEGY 

The LAE candidates that we aim to confirm in this pa- 
per are selected from our spectroscopic survey, which is 
presented fully in Paper II. In summary, the search uses 
the Inamo ri-Magellan Areal C amera and Spectrograph 
(IMACS; pressler et al.||2006[ ) on the Magellan-Baade 
Telescope, 'lb obtain blind, narrowband spectra we use 
a "Venetian blind" slit-mask (see Figure 2 of Paper I 
for a schematic) and a 134 A-wide (FWHM) blocking 
filter centered in the 8200A OH-free atmospheric win- 
dow. The slits on this mask cover an area of 55 arcmin^, 
corresponding to a volume of 1.5x10'* Mpc"^ for Lya at 
redshift 5.7. Emission lines that are detected are as faint 
as 2.5xl0~*® erg s~^ cm~^. Two fields were surveyed: 



one i n the center of the COSMOS field (Scoville et al. 
20071), and the other in the Las Campanas Infrared Sur- 



J, , . 

wytLCIRS) 15H field (|Marzke et al. 1999). In this pa- 
per we focus on the COSMOS held, as the 15H field has 
only limited spectroscopic followup at this time. A more 
complete luminosity function, with data from both fields 
combined, will be the subject of a future paper. 

We identified 105 LAE candidatefl!] in the IMACS 
search of the COSMOS field. These objects are seen as 
single emission lines with no detected continuum (or, in 
a few cases, a possible weak continuum detection only to 
the red of the line). An additional 130 objects with blue 
continuum or multiple emission line^were identified as 
foreground galaxies in the search data. 

Identifying the LAE candidates as genuine LAEs or 
foreground galaxies is not trivial. Since the positions of 
galaxies within the slits are not known, a large area on 
the sky must be searched for an optical counterpart. As a 
result, the normal approach of requiring an optical non- 
detection to distinguish foreground galaxies from LAEs 
is not applicable. This ambiguity is illustrated in Figure 
[l] which shows broadband and narrowband thumbnail 
images of our six confirmed LAEs (which we will discuss 
in q4). In most cases there are no objects falling within 
the slits at the spatial position of the LAEs, but in some 
cases there are objects nearby or outside the slits. Sim- 
ilarly, only 37 of the 105 LAE candidates have entries 
within a 2" radius in the COSMOS photometric redshift 
catalog]^ Of these, 20 have photometric redshifts that 

^ We have added one object that was dropped from the sample 
in Paper II, bringing the number of LAE candidates from 104 to 
105. 

* Multiple emission lines can be present in the 134 A of search 
spectrum when we detect Ha and [N II] A6583 A, or [O III] AA4959, 
5007 A. 

^ The photometric redshifts used in this paper are 
taken from the COSMOS Photometric Redshift Cat- 
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Fig. 1. — Thumbnail images centered on the locations of the confirmed LAEs presented in |4] Image are 5" on a side, and the dotted fines 
show the location of our slit. The Subaru (B, V, r, i, NB 816, and z), UKIRT (J), and CFHT (Ks) images are taken from the COSMOS 
image cutout service, hosted by the NASA/IPAC Infrared Science Archive. The i— dropout near 20.5+7-0.07 (see text) is circled in red. 



place an emission line in our 8200A bandpass and 17 do 
not. 

Clearly, the faintness of our LAE candidates means 
that spectroscopic foUowup is required to determine if 
the detected emission line is Lya, or an emission line 
from a foreground galaxy. Under this approach, we can 
detect other emission lines that identify the galaxies as 
foreground emitters, and also observe the "smoking gun" 
signature of Lya emission: an asymmetric line profile. In 
the sections below we describe our spectroscopic followup 
observations. The slits used to observe LAE candidates 
always have the same width, coordinates, and position 
angle (PA) as the search slits, regardless of potential op- 
tical counterparts. 

3. OBSERVATIONS 

3.1. IMACS Spectroscopic Followup 

Spectroscopic followup was carried out on March 23- 
26 2009, using the In amori Magellan Areal Camera 
(IMACS; [Dressier et~al. 2006) on the Magellan Baade 
Telescope! Conditions were clear, and the average see- 
ing over the three nights was 0.6". A slit width of 1.5" 
was used to match the 2008 search observations, and the 
200 line mm~^ grating provided a 2 A pixel" scale to 
go along with the 0.2" pixel"! g^g^jg q£ ^^le IMACS f/2 
camera. The full spectral range covered approximately 
4000 - 9000 A. 



alog (Fall 2008; version 1.5). which can be found at 
|http://irsa.ipac.caltech.edu/Missions/cosmos.html[ 



The data were reduced using the COSMOS (Carnegie 
Observatories System for Multi Object Spectroscopy) 
software package ^ These spectra allowed us to identify 
Ha, [O III], and~H/3 emission from galaxies. However, 
at this S/N and resolution neither the [O II] doublet nor 
the Lya line profile can be resolved. Therefore, both 
[O II] and Lya emitters appear as single emission lines 
in our IMACS followup spectroscopy. These single line 
objects are prioritized for higher-resolution spectroscopic 
followup with DEIMOS (see below). 

3.2. DEIMOS Spectroscopic Followup 

Medium resolution followup spectroscopy with 
DEIMOS ( jFaber et al.|[2003, ) on Keck II was carried out 
on March 7-8, 2010 and January 27-28, 2011. Further 
observations scheduled March 3-4 2011 were lost entirely 
to poor weather. In addition to the high-priority single 
line emitters from our IMACS followup, the slit-masks 
included LAE candidates which were not targeted for 
IMACS followup, and those that were targeted but 
not recovered. (Ultimately, it was not unusual for 
objects that we failed to recover in IMACS followup 
to be robustly detected under good conditions with 
DEIMOS.) The observing conditions and mask positions 
are listed in Table [T] In the 2010 observations, poor 
seeing impeded our ability to recover the faintest sources 
and to distinguish [O II] from Lya in low S/N emission 
line detections. Consequently, the optimal mask design 
in 2011 included substantial overlap with the previous 

http:/ /obs. carnegiescience.edu/Code/cosmos 
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TABLE 1 
Observation Summary 



Date 


Instrument 


Mask Name 


Mask RA 


Mask Dec 


Mask PA 


Exposure Time 


Seeing 


Conditions 














(hours) 


(") 




23-26 March 2009 


IMACS 


lOhconOg 


10:00:42.90 


02:11:00.0 


90 


16.0 


0.35-0.8 


Clear 


7 March 2010 


DEIMOS 


A 


10:00:29.19 


02:09:33.7 


95 


6.3 


0.9-1.6 


Clear 


8 March 2010 


DEIMOS 


B 


10:00:22.10 


02:01:38.0 


95 


6.0 


1.2-1.7 


Cirrus 


27 January 2011 


DEIMOS 


D 


10:00:22.97 


02:09:28.9 


85 


6.5 


0.6 


Clear 


28 January 2011 


DEIMOS 


F 


10:01:11.46 


02:10:27.4 


106 


6.3 


0.8 


Clear 



Note. — Coordinates and position angle (PA) are given for each observed mask. The slit PA, which differs from the mask PA given 
here, was 90 degrees for every object. Slit-widths were 1.5". 



year. Two LAEs identified on the March 2010 mask 
A (51.5-1-0.98 and 59.5+5-0.75) were also included on 
mask D in 2011. 

Observations were made with the 830G grating, us- 
ing slit-widths, position angles, and locations matched 
to the search data. The data were reduced using the 
DEEP2 DEIMOS Data PipelineF] with an updated op- 
tical model for the 830G grating (P. Capak, private com- 
munication). The January 2011 data were flux calibrated 
using observations of spectrophotometric standard stars, 
taken through a 1.5" slit at the parallactic angle. The 
stars u sed were G191B2B, GD50, Feige 66, F eige 67, and 
Hz 44 (|Okel|1990l IMas sey fc Gronwall|[T990| , with refer- 
ence data taken from the ESO spectrophotometric stan- 
dard star databasj^ Sensitivity functions derived from 
each of these stars agree to within 12%, and when Feige 
66 is excluded agreement is better than 3%. In addition 
to the default pipeline extraction, the January 2011 spec- 
tra are extracted in 1'.'8 apertures. This aperture size was 
chosen to match the spatial extraction in Paper II, and 
is appropriate given the good seeing in both the search 
data and January 2011 foUowup. No flux calibration was 
applied to the March 2010 data because variable seeing, 
at times larger than 2", and cirrus during some of the ob- 
servations, implied significant count-rate differences from 
frame to frame. 

Line fluxes were measured from the extracted January 
2011 spectra, and compared to the fluxes measured from 
the IMACS search data (and used in Paper II) . The root 
mean square deviation is approximately 50%, even for 
bright sources. This scatter likely results from the fact 
that the LAE candidates are not necessarily centered in 
the slits, so small astrometric errors or different seeing 
will induce varied slit-losses from one observation to the 
next. In ^ we take this scatter into account when we 
derive the luminosity function. For consistency, all line 
fluxes in this paper are taken from the search data. 

4. CONFIRMED LAES 

Six LAEs are identified from the DEIMOS observa- 
tions. Their coordinates, sizes, luminosities, velocity 
widths, and redshifts are listed in Table [2] In Figure 

twe show three cases of high quality Lya spectra, where 
e asymmetric line profile makes the identification un- 
ambiguous. Little discussion is needed for these three 
LAEs, although we do note that the two brightest LAEs 
51.5-1-0.98 and 59.5+5-0.75 are clearly spatially resolved 
on mask D, where they were observed under 0'.'6 seeing. 

http://astro.berkeley.edu/~cooper/deep/spec2d 

http: / / www.eso.org/ sci / observing/ tools /standards /spectra. html 



When an emission line is not clearly asymmetric, more 
care is needed to determine if it is Lya. Three such 
cases are identified, and are shown in Figure [3] and listed 
in Table [2j For these galaxies, other checks are pos- 
sible. First, we measure the S/N at the locations of 
[O III] AA4959, 5007 A, [O II]A3727 A, and H/3 for the 
expected foreground populations. In these cases, we find 
that II/3 and [O III] A4959 A are easily ruled out because 
other stronger emission lines would be detected. Like- 
wise, the [O II] doublet is ruled out by the line profiles 
shown in Figure [3] However, we can not rule out the 
case where we have detected Ha or [O III] A5007A, and 
all other lines are below the detection threshold in our 
DEIMOS data. While DEIMOS is efficient at 8OOOA, it 
is about a factor of two less sensitive at blue wavelengths 
where we would detect [O III] from an Ha emitter or 
[O II] from an [O III] emitter. Consequently, our spectra 
do not verify that the emission lines shown in Figure |3] 
are genuine LAEs. 

The lack of any imaged counterparts within the slits 
for 57.5+70.82 and 47.5+6-0.90 (Figures fl] and js]) sug- 
gests that emission lines may indeed by Lya. We next 
ask whether the implied high equivalent width would rule 
out Ha or [O III] emission. For an optical counterpart to 
be undetected in the Subaru i-band and the ACS /^m it 



must be fainter than about 27th magnitude (AB; |Henry 
et al.||2009[ [Taniguchi et al.||2009| . This gives observed 



frame equivalent width limits of > 160 and 210 A for 
these two candidates. If, on the other hand, the correct 
optical counterpart is one of the objects detected, but 
outside the slit, we would still expect a large equivalent 
width because the true line flux would be significantly 
attenuated by our slit. Nevertheless, using the COS- 
MOS NB816 catalogs which we derived in Paper II, we 
find that 18-23% of Ha and [O III] lines have equivalent 
widths higher than these thresholds. Since the extrap- 
olated counts of Ha and [O III] emitters from Paper II 
are substantial, the high-equivalent width subset of these 
foreground populations could exist in similar numbers to 
the LAEs. Nevertheless, it is unlikely that both 57.5+7- 
0.82 and 47.5+6-0.90 are foreground interlopers. 

One LAE candidate in Figure [3] has not been discussed 
so far: 20.5+7-0.02. Like the two other uncertain candi- 
dates, we can not rule out Ha or [O III] from the spec- 
trum alone. (Although, we note that this object was ob- 
served under very poor seeing in March 2010, so future 
observations should confirm or rule out this LAE candi- 
date.) However, the /814 postage stamp contains mul- 
tiple objects. We explore this further by examining the 
multi- wavelength photometry from the COSMOS survey. 
Thumbnail images are shown in Figure [l] remarkably, an 
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Fig. 2. — Spectra of high confidence, confirmed LAEs are shown in one dimension (top row) and two dimensions (middle row; dispersion 
in the horizontal direction). The bottom row shows 5" /gi4 cutouts from the COSMOS ACS observations (Koekcmoer ct al. 2007). The 
detection threshold in these im ages depends on the sizes of galaxies, but we can estimate that undetected sources have J814 ^ 27 ^Taniguchi| 
|et al.|2009[ [Henry et al.|2009^ . The dashed lines show the location of our 1'.' 5- wide slits. 



TABLE 2 
Robust and Tentative LAEs 



ID 



RA'' 
(J2000) 



Dec 
(J2000) 



Flux 



Luminosity Mask 



(10~^** erg s"-^ cm"^) (10"'"' erg s-^ ) 



FWRMi^ 
(") 



Redshift 



(km s-i) 



51.5-1- 


0.98 


10:00:49.766 


02:10:06.448 


9.97+ 1.34 


3.72 


A&D 


1.13 


5.730 


163 + 70 


40.5-8- 


0.55 


10:01:35.446 


02:07:22.911 


5.05 ± 0.63 


1.83 


F 


0.94 


5.666 


200 + 35 


57.5+7- 


-0.82 


09:59:59.604 


02:11:37.826 


4.07 + 0.58 


1.51 


D 


0.69 


5.717 


187+ 132 


59.5+5- 


-0.75 


10:00:11.686 


02:12:07.171 


10.0 + 0.59 


3.68 


A&D 


0.95 


5.701 


184 + 66 


47.5+6- 


-0.90 


10:00:07.145 


02:09:07.808 


5.41 + 0.70 


2.03 


D 


0.59 


5.749 


< 140<^ 


20.5+7- 


-0.07 


10:00:00.691 


02:02:24.636 


7.16 + 0.88 


2.67 


B 


1.22 


5.730 


< 140<i 



^ The position accuracy in the R.A. direction (along the slit) is 0.3" (RMS). 

^ Spatial FWHM are seeing convolved and measured in the DEIMOS foUowup spectra. (Accurate size measurements are not possible with 

the low S/N in the search data.) The seeing on each mask is given in Tabic [l] 

^ Velocity width of a truncated Gaussian model fit to the line profile, as described in 
These line profiles are consistent with being unresolved. 

object I'.'S east of our detected emission shows the ap- 
pearance of an dropout (circled in Figures [l]and|3]pj 
If the detected emission is uhiniately confirmed to oe 
Lya, we may have observed a fainter companion galaxy. 
However, the bluer foreground galaxy lying near the cen- 
ter of the thumbnail image is more closely aligned with 
the detected emission, but its redshift is uncertain. Due 
to its faintness, this galaxy is not included in the COS- 
MOS photometric redshift catalog. 

In summary, we find three certain LAEs with clear 
asymmetric profiles and three plausible LAEs, where the 



blue sensitivity of our foUowup spectroscopy is insuffi- 
cient to rule out all cases of Ha and [O III] A5007A. In 
the sections that follow we derive Lya LFs for two cases: 
the "3 LAE LF," where only the robustly identified LAEs 
are included, and also the "6 LAE LF" , where both the 
robust and plausible LAEs are included. 



LAE AND FOREGROUND GALAXY IDENTIFICATION 
RATES 



Although the i— dropout candidate in question was clearly 
inside our DEIMOS slit, no line emission is detected. This is not 
surprising, given the large UV-luminosity of this galaxy {Mjjy ^ 
—22.0 if it is at redshift 5.7), and the observation that Lya emission 
is not typical of such bright i— dropouts jStark et al.|2011[ l. Indeed, 
the NB816 photometry is consistent witn no line einission. 



In order to derive the Lya LF, we must address the 
question of how many LAEs remain in our survey that 
have not yet been spectroscopically confirmed. To do 
this, we supplement our IMACS and DEIMOS spec- 
troscopy with COSMOS photometric redshifts in order 
to identify as many foreground objects as possible. 

To use the photometric redshifts, we must first de- 
fine what we refer to as a foreground photometric match. 
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Fig. 3. — Same as Figure[2] but for less certain Lya identifications. The i— dropout near 20.5+7-0.07 (see text) is circled in red. 



As can be seen in Figure [T] nearby objects are present 
even around the robustly ioentified LAEs, so we can not 
identify foreground interlopers by the presence of opti- 
cal counterparts. However, the situation is significantly 
improved with t he COSMOS photometric redshifts (Il-l 
bert et al. 2009), because the galaxy templates contain 



emission lines, and the NB816 band is included in the 
photometric redshift derivation. The consequence of this 
approach, is that- if line emission is significant in the 
NB816 band- the photometric redshifts become "tuned" 
to the precise redshifts of our foreground populations. 
For the NB816 filter, these redshifts are: z = 1.17 - 1.2 
for [O II] A3727 A, z = 0.62 - 0.69 for [O III]AA4959, 
5007 A, and H^, and z = 0.20 - 0.25 for Ha and [S II] 
AA6716, 6731. To find photometric matches with all of 
our emission line galaxies, we search to a radius of 2" 
from the slit-center. This gives 20/105 matches among 
the LAE candidates, and 77/130 matches among objects 
identified as foreground galaxies in the search data. With 
the caveat that this technique cannot distinguish [O III] 
A5007 A from [O III] A4959 A or H/3, and possibly Ha 
from [S II] AA6716, 6731 A, spectroscopic foUowup shows 
that photometric matches provide the correct redshift 
identification in 24/24 observed cases 

In Table [3] we list the counts of foreground galaxies, 
LAEs, and unidentified objects by fiux bins. Line fluxes 
are defined in the first column, and are established to 
correspond to 0.2 dex bins centered at log{L/erg s"^) = 
42.1,42.3,42.5... for Lya at redshift 5.7. The total num- 
ber of emission line galaxies in the survey is given in 
column two, and broken into LAE candidates and fore- 

Many more foreground objects were followed up at the posi- 
tion of the photometric match and not the position of the IMACS 
search slit. Here, we only count the objects where the foUowup slit 
was at the same position as the search slit. 



ground galaxies (both from the search data) in columns 
three and four. Column five lists all foreground emis- 
sion lines that we have identified to date, including ob- 
jects that were initially LAE candidates but have fore- 
ground photometric matches or spectroscopic identifica- 
tions. Columns six and seven list the counts of confirmed 
LAEs for the six and three LAE cases, respectively. The 
final column (eight) lists the remaining unidentified LAE 
candidates. The sum of these unidentified object counts 
plus the six likely confirmed LAEs represent an absolute 
upper limit on the counts of faint LAEs in COSMOS. 
For accounting purposes, it should be noted that this 
conservative upper limit, plus the identified foreground 
galaxies column are equal to the iVobj given in the second 
column. Additionally, in these binned counts we have ex- 
cluded ten emitters fainter than log(F / erg s~^ cm~^) 
— —17.57 and two brighter than —15.57. 

We next aim to infer, from our spectroscopy, the 
number of LAEs and foreground galaxies that we ex- 
pect to find among the unidentified line emitters in Ta- 
ble [3j For this inference, we consider only the spec- 
troscopic identifications of objects without photometric 
matches, as the unidentified objects, by definition, ex- 
clude these matches. Additionally, we examine only the 
three faintest fiux bins given in Table [S) The brighter 
bins contain only small numbers of LAE candidates, and 
none of these objects that lack photometric matches have 
been observed with DEIMOS. Table |4 quantifies these 
foreground and LAE identifications. While the numbers 
of galaxies with spectra remain small, no specific popula- 
tion of foreground emitter or LAEs appears to dominate. 
(Although we will discuss the significant presence of H/3 
emitters in Appendix [B|) . 

In order to calculate the number of LAEs that are likely 
present among the unidentified objects, we must recog- 
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TABLE 3 
LAE AND Foreground Counts 



logF 




LAE candidates 


Search Data Foregrounds 


Identified Foregrounds 


LAEs (all) 


LAEs (secure) 


Unidentified 


-17.47 


42 


32 


10 


19 


1 





22 


-17.27 


54 


26 


28 


39 


2 


1 


13 


-17.07 


45 


26 


19 


33 


3 


2 


9 


-16.87 


30 


9 


21 


24 








6 


-16.67 


16 


3 


13 


14 








2 


-16.47 


14 





14 


14 











-16.27 


13 


1 


12 


13 











-16.07 


6 





6 


6 











-15.87 


2 





2 


2 











-15.67 


1 





1 


1 












Note. — Number counts are given for the various categories of objects in the search data, as well as after followup observations. Columns 
two tlirough four are based strictly on the search data, while tlie remaining columns talce our spectroscopic followup results into consideration. 
The final column gives the number counts of LAE candidates wliich have no spectroscopic or photometric identification. 

and upper limits (blue triangles). Under this approach, 
the lower limit is derived by including only the three ro- 
bustly identified objects. The conservative upper limit, 
on the other hand, also includes the three tentatively 
identified LAEs and all of the unidentified objects in Ta- 
ble |3] These discrete LFs are plotted in Figure |4j No 
slit-loss corrections are applied at this stage; they will be 




42.0 42.2 42.4 42.6 42.8 43.0 43.2 43.4 43.6 

Log(L / erg s"^ ) 

Fig. 4. — The discrete luminosity functions derived from our 
blind spectroscopic sample are shown in red (inferred from the 
low, 3 LAE and high or 6 LAE samples). The more conservative 
limits are shown in blue. For compari son, the brig ht-end data from 
the Suba ru Deep Field (SDF Shimasaku et al.'2006 Kashikawa et| 



[al 



2011| and COSMOS (Murayama et al. 20U7) arc shown, ihe 
f dat a shown here e xclude the tamtcst candidates reported by 



[Kashikawa et aL] | |2011[ l, where the spectroscopic confirmation rate 
is low. 

nize that these objects consist of two sub-classes. First, 
there are 27 galaxies where the low-resolution IMACS 
followup confirmed a single emission line, but was un- 
able to resolve the [O II] doublet. Followup spectroscopy 
with DEIMOS showed that these objects are predomi- 
nantly [O II] and Lya, but, to date, not all of them have 
been targeted for higher-resolution spectroscopy. On the 
other hand, the many galaxies which have no IMACS fol- 
lowup are observed to be a mix of Lya and all of the other 
bright rest-frame optical emission lines. These two sub- 
categories of emission line galaxies populate the number 
counts of unidentified emitters in Table [3j but the for- 
mer has a higher LAE fraction. Therefore, we split these 
counts, and give LAE fractions appropriate for each sub- 
category in Table |4j The total numbers of inferred LAEs 
are calculated for both the 3 LAE LF and also the 6 LAE 
LF. 

6. THE Lya LUMINOSITY FUNCTION 

6.1. Direct Calculation 

With the counts of LAEs given in Tables [3] and |4] a 
direct derivation of the LF is straightforward?^ The vol- 
ume subtended by the slits in our search mask is 1.5 x 10"' 
Mpc'^. Correcting these counts for the detection com- 
pleteness (AppendixKl), we derive the conservative lower 



addressed §6.2 

In order to derive a meaningful constraint on the 
faint-end slope of the LF, data that cover the bright- 
end and the knee are needed. For this calculation. 



we draw on LAE c atalogs available in COSMOS (Mu- 
rayama et"ar||2007[) a nd the Subaru Deep Field (SUF; 
Kashikawa et al.||2(Jll K. Shimasaku, private communi- 
cation). While spectroscopic observations are more sen- 
sitive to low equivalent width emission, no systematic 
correction is needed, because the contribution from these 
sources is negligible. No LAE candidates are seen with 
(convincingly) detected continuum in the IMACS search 
data, implying rest- frame equivalent widths, Wq ^7—10 
A. Therefore, the combination of these two different sur- 
vey methods is valid; 

The COSMOS ( [Murayama eFaT] |2007l ), and SDF 
(Kashikawa et al. 2U11) samples are derived in similar 
ways. Both are selected from Subaru/SuprimeCam ob- 
servations through the NB816 filter. However, the color 
cuts used to select the LAEs are not the same: Mu- 
rayama et al. require iz — NB816 > O.";^ while the 
SDF LAEs are selected with i-NB816 > 1.5. Compari- 
son of these catalogs shows only a small difference. All of 
the galaxies in the SDF catalog have iz — NB816 > 0.67, 
and only two have iz - NB816 < 0.70. Conversely, 83% 
of the Murayama et al. (COSMOS) catalog meets the 
SDF color cut. Since the suitability of one color cut over 
another is unknown, we make no corrections, but note 
that this level of systematic uncertainty is present in the 
data. 

For the SDF sample, we calculate the line fluxes and 
Lya luminosities, taking into account the non-uniform 
transmission of the NB816 filter (see Figure [5|. These 
quantities are calculated from the 2" aperture z— band 
and NB816 magnitude s, following the prescription in 



Kashikawa et al. 



(2011). When redshifts are unknown. 



The weighted average magnitude, iz, is commonly used to 
represent the continuum under the NB81 6 band, and is der ived 
from i and z fluxes: f j^ = 0.57/i + 0.43/z ( [Murayama et al.|2007[ 
Takahashi et al.|2007| . 
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TABLE 4 

LAE AND Foreground Rates among line emitters with no photometric matches 



log F Nspec No ID Lya (all) Lja (secure) Ho H/3 [O III] [O II] Other'' Inferred LAEs (all) Inferred LAEs (secure) 



Line Emitters with no IMACS foUowup 



-17.47 


8 


19 


1 





12 1 2 1 


2.1 





-17.27 


7 


9 


2 


1 


3 1 1 


2.3 


1.1 


-17.07 


6 


6 


3 


2 


1 3 1 


1.8 


1.2 


IMACS confirmed single line emitters 


-17.47 


2 


3 








2 








-17.27 


2 


4 


1 


1 


10 


2 


2 


-17.17 


3 


3 


2 


1 


1 


2 


1 


Total LAEs (Confirmed and Inferred) 


-17.47 












3.1 





-17.27 












6.3 


4.1 


-17.17 












6.8 


4.2 



The other foreground galaxies which we identify are two H7 emitters and a galaxy with blue continuum but an uncertain redshift. 



they are fixed to z — 5.70 so tliat the Une is assumed to 
fall at the center of the filter bandpass. A single aperture 
correction is applied to all galaxies in this sample, which 



we es timate by comparing SExtractor ( Bertin & Arnouts 
1996| ]V[AG_AUTO magnitudes to the 2" diameter mag- 



nitudes. This approach is preferable to simply adopting 
the MAG_AUTO magnitudes, which are unreliable at low 
S/N. From this comparison, we find a median aperture 
correction of -0.38 magnitudes to the NB816 photometry, 
indicating a +0.15 dex correction to the Lya luminosity. 

Considerable efforts to follow up the SDF sample have 
confirmed 46/8 9 LAE candidates arid identified four 
[O Ilj emitters ( Shimasaku et al^|2006 Kashikawa et al 



20111. The spectroscopic conhrmation is highly complete 
at the bright end, but falls off to fainter luminosities. 
Excluding the [O II] emitters, and dividing the sample 
at L(2") = 1042-6 erg s-\ 31/36 (86%) of the brighter 
LAE candidates are confirmed, but only 15/39 (38%) of 
the fainter candidates are confirmed. Since narrowband 
imaging is subject to the selection of spurious (no emis- 
sion line) sources at faint magnitudes, we restrict the 
sample to L{2") > lO^^'^ erg s~^ where the confirmation 
rate is high. The five unconfirmed LAE candidates with 
these bright luminosities are included, as their uncon- 
firmed status is likely the result of slit-mask design. 

We calculate the differential luminosity function for 
the SDF LAEs by determining the effective volume over 
which each object can be both detected with NB816 (2") 
< 26, and selected with i - NB816 > 1.5. Within 
the 725 arcmin^ of the SDF, these volumes range from 
1.7 X 10^ Mpc^ for L - lO''^-^ erg s'^ to 2.3 xlO^ Mpc^ 
for L > 10"*"^ erg s"^. For the five LAEs that are not 
spectroscopically confirmed we adopt an effective vol- 
ume defined by the redshift range that corresponds to 
the FWHM of the NB816 filter: 1.9 xlO^ Mpc^. Next, a 
completeness correction for each object, as a function of 
its NB816 magnitude (K. Shimasaku, private communi- 
cation), is applied. The differential luminosity function 
from the SDF is plotted in Figure li] 

Combined with the SDF, COSMOS is an excellent 
complement, covering nearly 10 t imes the area to a shal- 
lower depth. Over 1.86 degrees^ , Murayama et al. (20071 
report 111 LAEs. Spectroscopic redshifts for this sam- 




8000 8050 



8100 8150 .8200 
Wavelength (A) 



8250 8300 



Fig. 5. — The NB816 throughput (dashed line) is non-uniform, 
so the emission line flux inferred from NB816 magnitudes will, on 
average, be underestimated. For comparison, we show the much 
squarer transmission profile of the filter used in this work (solid 
line) . 

pie are unavailable, so like the unconfirmed SDF LAEs, 
we adopt a volume spanned by the FWHM of the filter 
bandpass: 1.9 x 10^ Mpc^. Likewise, line fluxes and lu- 
minosities are calculated in the same way that was done 
for the SDF LAEs. However, we note that, for the COS- 
MOS sample, the best estimate of the total luminosi- 
ties are taken from 3" diameter apertures. In Figure |4] 
we show the differential luminosity function for the 65 
LAEs brighter than 10'''^ erg a~^. Although the precise 
completeness of this sample is not known, comparison to 
the SDF LF shows that the COSMOS LAE catalog is 
complete above this conservative cut. As we will show 
below, the slightly lower normalization of the COSMOS 
Lya LF is actually a shift in luminosity of about 0.05 
dex, which results from the assumption that all of these 
emission lines are centered in the NB816 bandpass. 

6.2. Maximum Likelihood Derivation 

So far, in our calculation of the Lya LF shown in Fig- 
ure [4j we have assumed that the observed luminosity, L, 
is a good approximation to the true luminosity Lq . How- 
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ever, in some cases, this approximation will fail. In the 
blind spectroscopic data, a luminous galaxy can be de- 
tected as a faint object if it falls outside of the slit. Simi- 
larly, with narrowband imaging through the non-uniform 
transmission of the NB816 filter, luminous galaxies can 
be observed with faint luminosities if the line is trans- 
mitted far from the center of the bandpass. Not only are 
the luminosities altered by slit and filter transmission 
effects, but the volume from which a galaxy is drawn 
also becomes uncertain. Relative to fainter galaxies, in- 
trinsically bright line emission can be detected over a 
larger volume because these galaxies can be observed fur- 
ther from the slit- or bandpass- center. In short, blind- 
spectroscopic observations and narrowband imaging sur- 
veys are subject to the same problems: when precise 
slit-positions or redshifts are unknown, true luminosities 
are uncertain and volumes arc not well-defined. 

A maximum likelihood approac h to the luminosity 
functioii para meter estimation (e. g. Sandage, Tammann, 
& Yahil 1979[ STY) offers a natural solution to this prob- 
lem, as it allows the flexibility to model these effects of 
uncertain luminosity and volume. With a Monte Carlo 
simulation (described below) we can determine the prob- 
ability k{LQ — > L) that a galaxy with some true luminos- 
ity, Lq, will be observed with some other luminosity L 

' "2008). The luminos- 



( Gronwall et al. 



ity function tha 



2007 Reddy et al. 



unction of L is then 



we observe as a 
a convolution of the true luminosity function with these 
kernels, k{Lo L): 



(1) 



Here, 4>{Lq) represents the Schechter parameterization 
of the luminosity function: (j){x)dx = 0*a;"e~^da; , where 
X — L/L* . In the limit where there is no scatter between 
Lq and L,k{LQ — > L) becomes a delta function, and 
Equation n] reduces to (piL) = (/)(Lo)- 

This treatment can be applied directly to the STY 
maximum likelihood parameter estimation, where we 
choose a and L* that maximize the likelihood of observ- 
ing the set of N galaxies in our sample, with luminosities 
Lf. 

Including the convolution shown in Equation [l| the like- 
lihood function becomes: 



£(a,L*) = n 



/ fc(Lo Li)(j){Lo)dLo 
\JdLj dLok{Lo ^ L)(f>{Lo) 



(3) 



Since (j)* appears in the numerator and denominator, the 
normalization is set by the number of galaxies that are 
observed. To calculate the normalization, while account- 
ing for the covariance between a, <j)* , and L* , we calculate 
the number of galaxies, iV', predicted by each model: 



N' 



dL j dLQk{Lo 



L)0(io) X Volume. (4) 



The appropriate survey volume here is the one that is 
matched to the simulation used to derive k(LQ — > L). 
In the following calculations, the survey volume is not 
a function of luminosity, because galaxies are uniformly 



simulated out to some maximum distance from slit-center 
(or wavelength from away from the bandpass-center). 
Under this approach, k{Lo — )• L) reflects the incom- 
pleteness to fainter galaxies when their simulated po- 
sitions/wavelengths render them undetectable. Finally, 
we include the (Poisson) probability of observing the N 
galaxies in our survey, given the model prediction of N': 



C{a,L*,(l)*) = £{a,L*) * P{N\N'). 



(5) 



The joint constraint from the three different surveys that 
we use in this paper is simply the product of the likeli- 
hood functions, £(a,L*,0*), for each. 

6.3. Implementation of the maximum likelihood 
parameter estimation 

Having outlined the steps needed to derive the best- 
fitting Schechter parameters, we next implement this 
method for both our spectroscopic data and bright-end 
data taken from the COSMOS and the SDF narrowband 
imaging surveys. For each of these three surveys we use 
Monte Carlo simulations to calculate the convolution ker- 
nels k{Lo -> L). With these quantities in hand, we find 
the a, L* and (p* that maximize Equation [5j 

6.3.1. Blind Spectroscopic Data 

The quantities that determine fc(Lo ~^ L) for blind 
spectroscopy are: (1) the photometric scatter, and (2) 
the observed sizes of the Lya emission in our O'.'S see- 
ing search data. As mentioned in S|3j our DEIMOS ob- 
servations, when compared to the IMACS search data, 
indicate a photometric scatter with an RMS of 40-50%, 
regardless of line fiux. The sizes of LAEs, on the other 
hand, require a more careful consideration. Continuum 
observations show that the majority of LBGs, including 
those with Lyg emission, have compact sizes at z ^ 6 
( Bouwens et al.| | 2004{ i Ferguson et al. || 2004l [ Henr y et al. 
^010||Malhotra et al.p01l| [Gronwall et al |:!i01l[ ). Addi- 
tionally, Hubble Space Telescope imaging ot Lya emission 
from redshift 3-4 LAEs suppor ts the scenario w he re Lyg 
emission is similarly compact (Bond et al. 2010; Finkel-I 
stein et aL|2011) . However, two ot the securely idcntifieiT 
LAEs are clearly resolved in the DEIMOS foUowup ob- 
servations, so we consider the impact of size. 

Because our sample is too small at this time to adopt 
an empirically motivated size distribution, we consider 
three cases to illustrate the effect of LAE sizes. We 
assume constant sizes of O'.'S, 0'.'75, and I'.'O (observed 
FWHM) , and derive LAE parameter constraints (by the 
method described below). This test shows that, when 
larger sizes are adopted, the most likely values of (/)*, L*, 
and a are unchanged within the precision of our param- 
eter search (Alog(L*) = 0.01, Alog(0*) = 0.01, and Aa 
= 0.05). Only the likelihood contours are slightly altered 
when larger sizes are used. Relative to the constraints ob- 
tained when 0.5" FWHM are adopted, 1" sizes result in 
an extension of the likelihood contours to values brighter 
by 0.1 dex for log(L*), and steeper by 0.1 for a (68% con- 
fidence). For quantifying the uncertainties, we adopt the 
intermediate case where all LAEs have observed FWHM 
= 0.75" in the IMACS search data. 

With these assumptions, we determine fc(io i) by a 
Monte Carlo simulation where LAEs are randomly placed 
inside and near the slits, and noise is added. We consider 
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Fig. 6. — Left- The convolution kernels k{Lo — > L) describe the probability of a galaxy with luminosity Lq being observed v^ith some 
other luminosity L. While these probabilities are strongly peaked on the diagonal, where Lg ~ L, there is a non-trivial tail where galaxies 
with large Lq are detected with faint L. These are galaxies that fall partially outside of slits. This example is for the 6 LAE case, so the 
redshift identification incompleteness is apparent at L < 10*^'^ erg s~^. Right- A vertical slice of the left panel, taken for L = 10*^'^ erg 
s~^ shows the fraction of galaxies that are detected with this luminosity but in actuality have other luminosities, Lq. Again, the tail to 
high Lq comes from sources that fall partially outside the slit. 



an intrinsic luminosity range of Lq 



10 



41.5 



10*'^ erg 



The lowest simulated luminosity is fainter than the 
limit of our search, so this allows for photometric scatter 
into our selection. The high-luminosity cutoff is arbi- 
trary, but is justified since the volume densi ty of LAEs 



becom es vanishingiy small at this luminosity ( Murayama 
et al. [2007). This high— L cutoff determines the maxi 
mum distance from slit-center where LAEs are placed. 
A galaxy with Lq = lO^"* erg and a size of 0.75" 
(FWHM) can present detectable line flux even from a dis- 
tance of 1.7" from the slit-center. Therefore, for the blind 
spectroscopy, the volume used in Equation |4] is 3.4x10"* 
Mpc^. 

Since fc(Lo — ^ L) describes the probability of a galaxy 
with luminosity Lq being observed in our survey, it is 
appropriate to include survey completeness in the Monte 
Carlo simulation. This takes two forms. First, the de- 
tection incompleteness, which is derived in Appendix [A] 
is included. Additionally, since we are deriving LFs from 
inferred counts, we include a redshift identification in- 
completeness defined as the ratio of the confirmed LAEs 
(columns four and five) to the inferred LAEs (columns 
eleven and twelve) in Table |4j The bin with zero inferred 
LAEs is assigned a spectroscopic completeness of zero. 

Figure [6] shows k{LQ — ?> L) with the spectroscopic in- 
completeness included for the six-LAE case. The left 
panel shows that, in this two dimensional space, the 
probability is peaked on the diagonal because the most 
likely true luminosity, Lq, is the one that we observe, L. 
However, a tail of galaxies with high Lq and low observed 
L are also present when galaxies fall partly outside of a 
slit. On the left-hand side of this two dimensional space, 
at L < 10*^'^ erg s""'^, the spectroscopic incompleteness 
is apparent. The right panel of Figure [6 shows a ver- 
tical slice of the left panel to illustrate the fraction of 
galaxies at each Lq that are detected with L = lO'''^'^ erg 
s^*. Because of photometric scatter, there is a significant 



erg s ^ cm ^ for the 



probability of a galaxy having a fainter true luminosity, 
Lq, than is observed. It is worth noting that, if a set of 
Schechter parameters are assumed, we can use this slice 
in Equation [l] to get </)(L = 10*^-^ erg s~-^). 

6.3.2. Bright-end Narrowband Imaging Data 

To calculate k{LQ — >■ L) for the narrowband imaging 
samples, we must take two scenarios into consideration: 
when the redshifts of the emission lines are known, as 
they are for most (but not all) of the SDF LAEs, and 
when the redshifts are unknown, as is the case for the 
COSMOS narrowband imaging sample. When redshifts 
are unknown, emission lines with some Lq are placed 
randomly at wavelengths between 8050 A and 8250 A. 
Therefore, the appropriate volumes used in Equation |4] 
are 3.2 x 10^ Mpc^ for the SDF and 3.0 x 10^ Mpc^ for 
COSMOS. Additionally, a normally distributed flux error 
is added, with cr = 2 or 5 x 10~ 
SDF and COSMOS simulations. 

As with the blind spectroscopic data, these Monte 
Carlo simulations include detection completeness. While 
the COSMOS narrowband imaging sample is assumed to 
be complete at L > 10"''^ erg s~*, the SDF data must 
be corrected. Although the incompleteness for this sam- 
ple is a function of the NB816 magnitude, it varies slowly 
from 75 - 90% over the range that we are interested. Fur- 
thermore, the line luminosity forms a fairly tight correla- 
tion with the NB816 magnitude, with an RMS scatter of 
0.1 dex. Therefore, at each luminosity we adopt a com- 
pleteness appropriate for the mean NB816 magnitude of 
galaxies with that luminosity. 

When the redshifts of the emission lines are known, 
a different set of kernels is used in Equation |3] In this 
case, the simulation differs from the case of unknown red- 
shifts, because the observed L have been corrected for 
the transmission through the filter curve shown in Fig- 
ure [5j Therefore, in this simulation there is no scatter 
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Fig. 7. — Best-fit LFs to the inferred LAE counts (high and low), shown in red, are compared to other derivations of the reds hift 5.7 
lumin osity f unction. Curves are dashed at luminosities where they are extrapolated beyond the data in their various studies. For the |OuchI] 
|et al.| | |2008l l; |Kashikawa et al.||2011j l and jHu et al.| ( |2010| l studies we show the fits for a = —1.5, which is intermediate to the values favored 
by our data. 



between Lq and L from the unknown filter throughput. 
In the absence of noise, these kernels would be a set of 
(5-functions. However, as with the case of unknown red- 
shifts, noise and detection incompleteness are included. 
For the purposes of determining the normalization, we 
adopt the scenario used for unknown rcdshifts in Equa- 
tion H for aU SDF LAEs. 

6.4. Results: Luminosity Function Constraints 

With the various k{Lo — > L) in hand, we evaluate 
Equation [5] for each survey, over a grid of q, L* , and 
(f)* . The final likelihood function is given by the product 
of the individual likelihood functions, £{a, L* for 
each survey. The parameters that maximize this quan- 
tity are given in Table [5j for both the 3-LAE and 6-LAE 
LFs. 

The best-fitting Schechter-functions are plotted in Fig- 
ure [7] One feature of maximum likelihood parameter es- 
timation is that the most likely model does not always 
run through the data points. This is indeed the case in 
Figure [7j The discrete points describing the COSMOS 
LAEs are derived by assuming that all objects are ob- 
served through the peak of the NB816 throughput curve. 
Consequently, on average, the true luminosities, Lq, of 
these LAEs will be brighter than what we have observed. 
Since we model these effects with our maximum likeli- 
hood estimation, we expect the Schechter functions plot- 
ted in Figure l7| to fall above the discrete data points for 
the COSMOSnarrowband imaging data. On the other 



hand, since we have taken the bandpass into account for 
the 85% of the SDF LAEs that are spectroscopically con- 
firmed, we expect the most likely Schechter function to 
closely follow the discrete SDF points in Figure [7j In 
fact, if we exclude the SDF data from our fit, the bright- 
end of the best-fit Schechter function falls close to the 
same location, approximately 0.05 dex brighter than the 
COSMOS narrowband imaging points shown in Figures 

tand [7] Correcting these points in Figures |4] and |7] for 
is average attenuation, would bring them in line with 
the SDF data and the Schechter function fit. 

Figure [s] shows the 68% and 95% confidence intervals 
for these parameter constraints. Clearly, owing to our 
small sample of faint LAEs, a large range of acceptable 
faint-end slopes can describe our data. However, as our 
spectroscopic foUowup becomes more complete, the lu- 
minosity function parameters will be more tightly con- 
strained. 

7. IMPROVING THE FAINT-END SLOPE CONSTRAINT 
WITH MORE FOLLOWUP SPECTROSCOPY 

It is not surprising that the luminosity function pa- 
rameter constraints which we show in Figure [8] allow a 
wide range of Schechter parameters. At this time the 
faint-end of our LF is based on only a few confirmed 
LAEs. Our continued efforts at followup spectroscopy 
will narrow the allowed parameter space by increasing 
the sample size. Furthermore, the addition of our 15H 
field, when spectroscopy is complete, will approximately 
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Fig. 8. — Likelihood contours show the Schechter parameter constraints from the MultisUt Narrowband Spectroscopic sample (this paper), 
combined with the bright-end narrowband imaging data from COSMOS and the SDF. The 68 and 95% confidence intervals are represented 
by the shaded contours for the 3 LAE LF and by the open curves for the 6 LAE LF. 



TABLE 5 

Most Likely Schechter Function Parameters 



Sample 


log(0*/Mpc~^) 


log{-^-*/crgs ^) 


a 


3 LAE 
6 LAE 


'-'•^^-0.52 
o oc+O.Sl 
'^•'-'^-0.56 


42.91«-i? 
42.96«:?^ 


_1 45+0.92 

-1 70+0'^3 
"-'-0.57 



double the number of LAEs over what is present in the 
COSMOS field. 

It is therefore of interest to address how strict the con- 
straints from our survey may ultimately prove. To ex- 
plore this, we create a mock catalog of galaxies with (/>* = 
lQ-3.35 Mpc"^ L* = 1042-96 ergs^i, and a = -1.70 
(our 6 LAE LF in Table [s] and the contour curves in 
Figure [s]) . These galaxies are distributed throughout the 
volume of both our COSMOS and 15H IMACS fields, out 
to 1.7" fro m sl it-center, similar to the Monte-Carlo simu- 
lations in §6.2[ Slit-losses, noise, and detection complete- 
ness are included. Finally, we assume a spectroscopic 
completeness of 90% at all luminosities, since achieving 
100% completeness is usually impractical. The resulting 
sample of LAEs in both fields is approximately 40 galax- 
ies. After regenerating the convolution kernels to include 
the appropriate spectroscopic completeness, we calcu- 
late the likelihood function, C{a,L*, cj)*) over the grid of 
a, L* , (f>* (leaving the bright-end data unchanged). Con- 
tours shown in Figure [9] illustrate the improvements in 
the constraints that we can expect to find from complet- 
ing our spectroscopic followup. For the increased sample 
size in this simulation, the 68% confidence interval on 
the faint-end slope is reduced by 50%. Likewise, under 
our current constraints, the space density of LAEs at 
L — 10^^ erg s~^ spans a factor of 10 (68% confidence), 
but should be decreased to less than a factor of three 
when both fields are 90% complete. 

8. DISCUSSION 

8.1. Comparison to previous measurements of the 
redshift 5. 7 Lya LF 

In Figure [7] we compare the LFs derived in this work 
to those presented in previous studies at redshift 5.7. A 
number of items are readily evident. First, the shape 



and L*) of the LF that we derive from th e SDF 



data is diffe rent from the shape obtained by .Kashikawa 



et al. (2011), from the same bright-end data. This varia- 



tion IS not surprising, given our different analysis. Most 
notably, we have adopted luminosities calculated from 
the NB816 and z-band magnitudes, whereas Kashikawa 
et al. use luminosities derived from their followup spec- 
troscopy. Given the uncertainties associated with esti- 
mating slit-losses, imaging-based measurements may be 
a better approximation to the total Lya luminosity- pro- 
vided that aperture magnitudes are corrected to total 
magnitudes, and the redshifts are used to calculate the 
true luminosity observed through the NB816 filter. 

Figure [7] also compares our bright-end LF deriva- 
tion with the similar redshift 5.7 Lya LF from the 



Su baru/X MM-Newton Deep Survey (SXDS; Ouchiet 



20081). While the agreement between the plotted 
lechter functions is good, this may be fortuitous. Lim- 
ited spectroscopic followup in the SXDS could allow sig- 
nificant contamination, since narrowband imaging can 
sometimes select LAE candidates which have no line 
emission at all. This point is illustrated clearly in the 
color magnit ude plots used to select LAE candidates 
(Figur e 7 inJOuchi et a l.||2"008[ Figure 1 in Murayama 
et al.||2007l Figure 7 Tn Sliimasaku et al.||2006p. When 



the limit or the NB816 photometry is approached, there 
is significant scatter about the threshold that defines an 
emission line source. While the use of optical "veto" 
bands will eliminate some spurious emission line candi- 
dates, even this method could break down for faint galax- 
ies, where modestly red galaxy SEDs may be undetected 
at bluer wavelength s. This contamination m ay account 
for the LAEs which [Kashikawa et al. (2011) do not de- 
tect in followup spectrosco py. Likewise, the fa intest con- 



firmed LAEs in the SXDS fOuchi et aL]|2008[ ) are about 
one magnitude brighter in MB816 tlian the limit of that 
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Fig. 9. — Likelihood contours show the Schechter parameter constraints that we can achieve when the spectroscopic foUowup is 90% 
complete in both of our survey fields (curves). An improvement relative to our current constraints (shaded) is expected. For simplicity, we 
compare only the LF derived from 6 LAEs to a simulated sample derived from the 6 LAE LF Schechter parameters. The bright end data 
arc unchanged from what we have used in this study. Contours are shown for 68 and 95% confidence. 



sample. More spectroscopic followup will be needed to 
deterraine the reliability of the fainter Lya candidates 
selected from narrowband imaging. Such observations 
would provide better constraints on the space density of 
LAEs at luminosities around the knee of the LF, and 
reduce the covariances shown in Figures [8 and |9] 

Two other measurements deviate from the LF that we 
derive here: 



Hu et al. 



(|2010|), and Paper L These bright 
and intermediate luminosity data would require a discon- 
tinuity from our new faint-end results. A Schechter func- 
tion can not be fit which simultaneously goes through 
these data and our new data. Therefore, it is clear that 
both the Paper I and Hu et al. samples are incomplete. 
In deed, incompleteness in t he Hu et al. result is discussed 



Kashikawa et al. 



(2011). Because Hu et al. require 
spectroscopic confirmation, but use shorter observations 
than Kashikawa et al., it is argued that the incomplete- 
ness is simply a failure to re-detect some real LAEs. On 
the other hand, for our pilot survey, comparison of our 
current completeness to the detection completeness de- 
rived in Paper I suggests that the prior estimate was 
overly optimistic. For the current data, near 100% com- 
pleteness is reached for detections of better than lOcr sig- 
nificance. Shifting the completeness curve for the present 
data (Appendix Kl) by a factor of four to five in fiux al- 
lows us to revise the completeness estimate used for the 
shallower data in Paper I. This exercise suggests that 
our previous LF was 30-70% complete, rather than near 
100%, as was assumed. Revising these completeness cor- 
rections would bring the Pap er I points more clo sely in 
line with our new LF and the Ouchi et al. (20081 result. 

In Paper II, our analysis of the foreground populations 
suggested that the remaining counts of LAEs followed a 
similar faint-end slope, with a ~ —2. With our followup 
spectroscopy, we now find that the most likely slope is 
slightly shallower, although a = —2 is well within the la 
measurement uncertainty. Nevertheless, we identify the 
source of this difference (for the COSMOS field) in Ap- 
pendix IB] Briefiy, the foreground emission line galaxies 
that we nave identified so far in our followup add up to 
the counts predicted by the foreground model adopted 
in Paper II. Since we expect that a significant fraction of 
the unidentified emission lines listed in Table |3] are also 



foreground objects, we expect these counts to be higher 
than the most likely model adopted in Paper II. Because 
of this, the fraction of LAEs among our detected emis- 
sion lines will be smaller than we previously estimated 
in Paper II. As we show in Appendix [B] the most likely 
cause of this underestimation is a strong presence of H/3 
emitters from an overdensity at z ^ 0.68 in the COSMOS 
field. 

Finally, we note that the cosmic variance for our faint 
LAE sample, at this time, makes a negligible contribu- 
tion to the error budget because the Poisson errors on 
our small sample are much larger than the cosmic vari- 
ance. On the other hand, at the bright end, the excel- 
lent agreement between the narrowband imaging samples 
from COSMOS, the SDF, and the SXDS is better than 
the approximately 20- 30% cosmic variance th at we ex- 
pect for these surveys ( Trenti fc Stiavelli]|2009 ) . 



8.2. Comparison to faint Lya LFs at lower redshift 

So far, we have compared our Lya LF to previously 
measured Lya LFs at redshift 5.7. By necessity, this dis- 
cussion focussed on the bright-end data, since no faint- 
end measurements existed prior to the one which we 
present here. In this section, we turn to lower-redshift 
spectroscopic surveys which measure LA Es to very faint 
luminosities at z 3 — 5. In Figure [TOl we show three 
measurements of the cumulative Lya luminosity func- 
tion. Each of these searches use bro adband spectr oscopy 
to detect Lya, either with lensing ( Santos et al. [2004 j 



through serendipitous detections ( Sawicki et al 
or by exce ptionally long expo sures with the Very 
TelescopeJIlauch et al.l 20081) 



Large 



Figure[TO|shows that the space density of LAEs spans a 
factor of several between these different measurements. 
Large statistical errors, coupled with differing analyses 
and a wide range of redshifts make it difficult to draw 
any conclusio ns. The highes t space density of LAEs is 
found by Rauch et al. (2008), who also have the great- 
est sensitivity to low surface brightness emission. This 
could mean that they measure brighter luminosities by 
including more low surface brightness flux, and also that 
they are able to detect extended, luminous LAEs that 
other surveys would miss. However, where Rauch et al. 
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Fig. 10. — Spectroscopic surveys measure the cumulative Lya luminosity functions to very faint luminosities at z > 3. These include a 
search for strongly lensed LAEs (Santos et al. 2004), a serendipitous survey (Sawicki et al. 2008), and a very long exposure with the Very 
Large Telescope ( |Rauch et al.|20U8i ). i'he shaded area show the 68 and 95% confidence intervals from our LF fits, where, for simplicity we 
have taken the lower bounds trom the 3 LAE LF and the upper bounds from the 6 LAE LF. 



overlap with our survey, the agreer aent is good. Lower 
den sities are found at z ^ 4 — 5 by Santos et aL (2004) 
and [Sawicki et aL| (|2008| ). At L = 10^^" erg s~' ' where 
they o verlap with the wider, shallower search from |Ouchi 
et al. (j2008), the density of LAEs found by Santos et al. 
and Sawicki et al. is several times lower. The origin of 
this discrepancy is unclear. In short, statistical uncer- 
tainties make it difhcult to determine if the faint end of 
the Lya LF evolves at z ~ 3 — 6. Larger spectroscopic 
samples of faint LAEs are needed. 

8.3. Comparison to UV-selected i— dropout samples 

Surveys for LAEs have have often been motivated by 
the idea that line emission can be more easily detected 
than continuum when galaxies become very faint. This 
rationale is especially true for our spectroscopic search, 
which includes lines to fainter luminosities than most 
other Lya searches. For example, at the median LAE 
luminosity in Table[2](i — 10^^'^ erg s""'^), a galaxy with 
rest Lya equivalent widtlj^ Wq, larger than 190 A will 
have MiJY > —18. At 2 ^ 6, this luminosity is beyond 
the limit of the z— dropout samples used to measure the 

^® Here, and in the discussion that follows, we assume that fx oc 
with /3 = -2.0. While some evidence exists fo r bluer slopes 
at the highest redshifts ("Bouwens et al.|200 9' '2010b). these claims 
are disputed (TKinlop et al. 2011 see also, IJ'inkelstein ct al. 2010). 
Until more conclusive evidence for bluer slopes is found, we adopt 
a flat UV-slope when converting from rest-frame 1350 Ato 1216 A. 



UV-luminosty function ( Bouwens et al.|[2007 ). Likewise, 
at the L = 10^^ *^ erg s~^ limit of our survey, a galaxy 
would need to have a rest equivalent width smaller than 
about 80 A to have a continuum flux-density which is 
bright enough for inclusion in the i— dropout sample. 

Of course, assessing whether LAE surveys are finding a 
significant population of previously unobserved galaxies 
is not straightforward. The Lya equivalent width of star- 
forming galaxies displays a wide range of values, with a 
tail t o quantiti es larger than 100 — 200 A (Shaplcy et al. 
20031 [ Shimasa kuet al.|[2006{ [Gronwall et ai.[|2007| [S'tax^ 
et al. [2010) |20lT| ). Consequently, a faint lAe could be 
the result oi a low equivalent width emission from a UV- 
bright galaxy, or conversely, high equivalent width emis- 
sion from an extremely faint galaxy. No obvious trends 
exist between Lya luminosity and equivalent width (see 



Henry et al. 2010), although foUowup of UV-continuum 
selected i— dropouts shows that fainter continuum is cor- 
related with stronger an d more common Lya emission 



(Stark et al. 2010 20111. Ultimately, since Lya equiva- 
lent width exhibits such a wide range of values, a statisti- 
cal comparison with the ensemble of i— dropout galaxies 
is most relevant. 

To determine whether the LAEs in our sample are 
likely to be detectable as faint i~ dropouts, we gener- 
ate a mock catalog of galaxies. Beginning wit h the UV- 
lumin osity function of z ~ 6 i— dropouts (Bouwens et al 
X 10-3 Mpc-3, M^y = -20.24, and a = 



2007 



= 1.4 
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— 1.74), we generate a sample of galaxies in a fixed vol- 
ume, with UV- luminosities extending to Mjjv — —14.0 
(beyond the depth required for all LAEs in our survey to 
be detected in the continuum). A rest- frame Lya equiv- 
alent width (Wq) is assigned to each galaxy in this vol- 
ume, drawing from an exponential distribution with a 
maximum of Wq = 2000 A. Motivated by the observa- 
tions shown in Stark et al. (2011; Figure 2), we adopt e- 
folding, or "characteristic" equivalent widths of 50 A for 
Muv > M'ljy, and 20 A for Muv < M'^jy. These scale 
factors approximately reproduce the distribution of rest 
equivalent widths and the fraction of sources with equiv- 
alent width less than 25 A. Since narrowband imaging 
surveys cannot detect LAEs with < 25 A, galaxies 
which are given these low equivalent widths under our 
exponential parameterization are excluded from the Lya 
LF. (Although we note that the Lya LF is not strongly 
dependent upon this cutoff.) With each galaxy assigned 
an Muv and Wq, we calculate the Lya luminosity of each 
i— dropout. To understand the impact of the current 
observational limit, we repeat this calculation including 
only galaxies brighter than Muv = —18. The resulting 
simulated LFs are plotted in Figure [TT] where they are 
compared to our measured constraints on the Lya LF. 

The simulated LFs shown in Figure [TT] illustrate the 
importance of surveys for faint LAEs. With the equiva- 
lent width distributions chosen above, the Lya LF sim- 
ulated with only galaxies having M(jv < — 18 (dashed 
curve) shows a significant deficit of LAEs relative to when 
the complete range of luminosities are included (solid 
curve). At the L — 10*^ erg s~^ limit of our survey, 
approximately 50% of LAEs should have Muy > —18. 
Galaxies as faint as Muy ~ —16.5 are needed to bring 
these two simulated LFs into agreement. This result 
clearly indicates that we have uncovered a population of 
galaxies whose continua are too faint to have been ob- 
served in the deepest space-based dropout surveys. On the 
other hand, the narrowband imaging surveys for brighter 
LAEs are not unco vering a previously unobserved pop- 
ulation. Figure 11 shows that, brighter than L = lO''^'^ 
erg s~^, the Lya luminosity function is dominated by 
galaxies with Muv < —18.0. 

Figure 11 shows that the simulated Lya luminosity 
function rails short of observations at the bright end. 
A scale factor of 30 A rather than 20 A for the bright 
{Muv < ^uv) *~ dropouts would bring the simulated 
luminosity function within the 68% confidence interval 
shown in Figure [TT] A larger scale factor is not implau- 
sible, as the sa mple of bright z— dropouts targeted for 
spectroscopy by [Stark et al. (2011 1 is not large enough to 
derive very precise constraints on the luminosity depen- 
dence of the Lya equivalent width distribution. Addi- 
tionally, it should be n oted that the equiv alent width dis- 
tributions measured by Stark et al. ( 2011 1 are for galaxies 
that are selected as i— dropouts. I he selection function 
for thes e galaxies is known to depend on Lya equ ivalent 
widths ( Bouwens et al.|2006 Stanway et al.|2007 ), so the 
"raw" equivalent width distributions presented by Stark 
et al. may not strictly represent the parent population 
of star-forming galaxies. 

One way to interpret the (incomplete) simulated Lya 
luminosity function of galaxies with Muv < 18 is as a 
prediction of the Lya LF that one could measure from 
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Fig. 11. — Simulated Lya LFs (solid and dashed curves) are 
compared to our LF constraints. The shaded regions show our 
68 and 95% confidence limits. For simplicity, the lower bounds 
are taken from the 3 LAE LF and the upper bounds are taken 
from the 6 LAE LF. The simulated LFs are derived by including 
galaxies to very faint limits {Mjjy = —14; solid), and also by lim- 
iting the simulation to luminosities that have been observed in the 
UDF and UDF-Parallels [Muv < -18; dashed). The latter case 
turns over, representing an incompleteness to faint LAEs among 
the i— dropout samples. At the limit of our MNS survey (10'*^ erg 
s~^), 50% of LAEs have UV luminosities fainter than the limit of 
current i— dropout samples. 



spectros copic foUowup of i— dropout galaxies (Jiang et 
al. 12011 ) . However, correcting this LF to the "true" Lya 
LF ot all galaxies is not possible. We could speculate 
that the characteristic equivalent widths may be larger 
for galaxies with Muv > — 18 (following the trend of 
more prevalent Lya emission among galaxies with fainter 
UV continuum luminosities). Under this scenario, the 
discrepancy between the simulated LFs shown in Figure 
[TT] would be amplified, because the faint-end slope of 
me "true" LF (simulated to a realistic low-luminosity 
cutoff) would be even steeper. Since we cannot know 
the equivalent width distribution of the faint, continuum- 
undetected galaxies, only dedicated emission line surveys 
will be able to measure the faint-end slope of the Lya LF. 

8.4. Lya line profiles 

The transmission of Lya photons is sensitive to the 
ionization state of the IGM, so the shape of the Lya 
emission line profile may offer a measure of the changing 
neutral hydrogen fraction. In a more neutral IGM, the 
Gunn-Peterson Lya damping-wing will lower the peak 
ffux of the line, effectively broa dening the emission and 
remov ing its asymmetric shape ( Ouchi et al.|2010 Dayal 
eraD[2008). 

Observational efforts have begun to search for a dif- 
ference in Lya profiles observed at redshifts 5.7 and 6. 5 



Ouchi et al. 2010 Hu et al. 2010 Kashikawa et al. 2011 1 



but the average line width is not significantly ditterent 
between the two redshift bins. However, Hu et al. show 
that the Lya line width disp lays a positive correlation 
with luminosity. In Figure [12] we combine our faint LAEs 
with the Hu et al. redshift 5.7 sample to extend the 
correlation to lower luminosities. As expected from an 
extrapolation of the brighter sources, our newly mea- 
sured line widths are among the narrower values. The 
significance of the relation is improved: under the Spear- 
man rank correlation test, the probability of the null hy- 
pothesis (no correlation) is decreased from 1.6 x lO^'^ to 
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Fig. 12. — The hya line width shows a positive correlation with 
luminosity. Here, cr is determined in the same way for all LAEs: a 
truncated Gaussian is convolved with the instrumental resolution 
and fit to the observed line profile. Two objects have line profiles 
that are consistent with being unresolved, and for these we plot 
upper limits that correspond to the instrumental resolution. 

2.6 X 10^^ when our six faint LAEs are included. These 
probabiUties, interpreted in the context of a Gaussian 
probability density function, increase the significance of 
the correlation from Scr to better than 3.5a. 

The luminosity dependence of the Lya line-width has 
important implications for how the effects of reionization 
on the line profile are interpreted. Comparisons made at 
different redshifts should also be made at fixed luminos- 
ity. This luminosity dependence has not been accounted 
for when entir e samples of redshi ft 5.7 an d 6.5 LAEs hav e 
been stacked ( Ouchi et al.||2010 Kashika wa et al.||2011 ). 
In a flux-limited sample, the redshift 6.5 LAEs will be, 
on average, more luminous. This could explain the ob- 
servation that the redshift 6.5 line profiles are marginally 
broader than those at redshift 5.7. 

The origins of the observed lumino sity dependence are 



unclear. One explanation offered by Hu et al. (2010) is 



that the more luminous LAEs may simply be more mas 
sive, with intrinsically broader Lya lines. Another pos- 
sibility lies in the kinematics and conditions of the cir- 
cumgalactic medium. While typical star-forming galax- 
ies have redsh ifted Lya emission and blueshifted P- Cygni 
absorption ( ,Shapley et al.|2063 Steidel et al.|2010 l, rare 
cases are seen wh ere blueshifted material emits rathe r 
than absorbs Lya ( |Erb et al.|2010[|McLinden et al.|2011[ ). 
In these cases, the Lya emission is both luminous, and 
broader than average. While Gunn-Peterson absorption 
will destroy blueshifted emission at redshift 5.7, the rela- 
tive strength of the emission component at zero velocity 
may still induce a correlation between line width and 
luminosity. Measurements of LAE masses and systemic 
redshifts can distinguish these scenarios, but for reioniza- 
tion era galaxies, the required observations will remain 
challenging until the James Webb Space Telescope oper- 
ates. 



8.5. The contribution of faint LAEs to the 
hydrogen-ionizing background 

At z ~ 6 the reionization of the intergalactic medium 
was complete, and the neutral hydrogen fraction in the 
IGM was negligible. Star-forming galaxies are usually 
assumed to produce the necessary ionizing background 
to prevent hydrogen recombination. However, observa- 
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Fig. 13. — Including galaxies fainter than those that we have 
observed increases the Lya luminosity density, making it easier 
for galaxies to maintain an ionized IGM. Shaded regions show the 
68 and 95% confidence limits on the Lya luminosity density, as 
a function of the low-luminosity cutoff. To simplify, the upper 
bounds are taken from 6 LAE constraints and the lower bounds 
are taken from the 3 LAE constraints. The critical Lya luminos- 
ity density required to maintain an ionized IGM at redshift 5.7 is 
Lp{Lya) = C, y. 3.0 x 10'"' erg Mpc~^. Horizontal lines mark 
C = 0.1,0.5,and 1.0. 



tional constraints at redshifts 6 and higher have yet to 
demonstrate that galaxies produce the requisite ionizing 



background (Bunker et al. 



2007[ 

too taint to' be observed are invoked 
2004; Salvaterra et 



2004 2010 Bouwens et al. 



2010a Uesch et al.|20l0| ) . Often, galaxies which are 



al|2011| 

s from b: 



lYan fc Windhorst 

_^ Finlator eral.|201ip , extrap- 
olating taint counts from brighter populations. In this 
paper, we address this shortcoming by detecting some 
objects from this hitherto undetected, postulated faint 
population ( as d emonstrated in [8.3). 

In Figure [13] we show the cumulative Lya luminos- 
ity density (as a function of the low-luminosity cutoff) 
allowed by our newly derived constraints. In order to 
determine whether LAEs can maintain an ionized IGM, 
we must infer an ionizing photon rate from the observed 
Lya photons, and then address whether this ionization 
rate can balance the recombination rate in the IGM. 
This threshold depends on the escape of both Lya and 
hydrogen-ionizing (LyC) photon s (/ldC and fLva\ Siana 



et al.||2010[ [Nestor et al1|2011 



Hayes et al. 201T), as 



well as the dumpiness of ionized hydrogen in the iGM 
(C = (j^jj jj)/('^H ll)^)- -^^ Paper I and Paper II we 
renormalized these uncertain parameters to their best 
guess values and grouped them into a single parame- 
ter, Q. In this way, C = 0-9 when C = 6, fLya = 0.5, 
and fhyC = 0-1 7 and the critical Lya luminosity den- 
sity required to maintain an ionized IGM at redshift 5.7 
is Lp{Lya) = C x 3.0 x 10-*° erg s^^ Mpc"^ Thresh- 
olds with C = 0.1, 0.5, and 1.0 are shown in Figure [T3[ 
This plot illustrates that, if C, is among the smaller values 
of 0.1 to 0.5, it may be possible to maintain an ionized 
IGM with only the galaxies that we have observed. How- 
ever, if C, is larger, even fainter galaxies are needed and 
the faint-end slope must fall among the steeper values 
allowed by our data. As we pointed out in Paper II, 
when faint-end slopes are shallow (a > —1-3), the Lya 
luminosity density grows very slowly when increasingly 
lower luminosities are included. Therefore, in the case 
of a shallow slope, it may be difficult for star-forming 
galaxies to reionize the IGM. 
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8.6. Properties of galactic building blocks at the end of 
reionization 

In Paper II we compared our density of faint LAEs to 
that of local galaxies, estimating that we had found three 
or four sites of active star-formation per L* galaxy today. 
The implication of this result was that we had identified 
the building blocks of Milky- Way type galaxies. In this 
work, we have refined our faint-end slope measurements, 
and shown that the most-likely slope is shallower than 
the one we adopted in Paper II. As a result, the space 
density of LAEs that we have identified is about a fac- 
tor of two smaller- 2 x 10~^ Mpc~^. However, not all 
galaxies have Lya emission. By comparison, the space 
density of observed i— dropouts is about 2-3 times higher 
than the LAE density that we observe. Furthermore, 
if z— dropouts were detected as faint as Muv = —16.5 
(our estimate of our faintest continuum luminosities from 
^8.3), their space density would be several times higher 
than the LAE density presented here. Since our LAEs 
belong to this class of postulated ultra-faint i— dropouts, 
we can infer that we have observed the building blocks 
of present day L* galaxies. 

While it remains difficult to determine the physical 
properties of LAEs from their Lya emission line alone, in 
Paper II we placed their properties in the larger context 
by comparison to other studies. For example, the angular 
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correlation o f the brighter LAEs presented in i 
( |2008j|2010| implies halo masses of M = 10^^ 
at redshitt 5.7. Presumably, the lower luminosity objects 
that we have found reside in less massive halos. Addi- 
tionally, in Paper I we estimated SFRs as low as 1 to 
a few Mq yr"""^ (if half of the Lya photons escape the 
galaxy) . If this star formation rate persists for a dynam- 
ical time scale of at least 50 Myr, then approximately 
10^ Mq of stars will have formed. These masses, SFRs, 
and redshifts are well matched to the predicted proper- 
ties of the galaxies that enriched the IGM w ith metals 
( [Oppenheinier et al.|[2009{ jMartin et al.||2010[ ). 

As noted m some ot the LAEs have spatially re- 
solved emission in our DEIMOS spectroscopy. This is 
not strictly surprising, as stacking of narrowband Lya 
images at 2: ~ 2 — 3 has revealed that diffuse, extended 
Lya e mission is a generic component LA Es and LBGs 



alike (Stcidcl et al. 2011). Additionally, Ranch et al. 
pOOS) find that, when exceptionally faint surface bright- 
ness limits are reached, a considerable fraction (40-70%) 
of redshift 3 LAEs are extended. Consequently, the ob- 
servatio n that LAEs are compact in narrowb and HST 
images (Bond et al. 2010 Finkelstein et al. 2011) does not 
imply that the Lya emission region is small, but rather 
that the sensitivity to low surface brightness emission 
is insufficient in the space-based images. Further stud- 
ies of extended Lya emission in larger samples may help 
quantify the evolution of the circumgalactic medium at 
high-redshift. 

9. SUMMARY AND CONCLUSIONS 

We have measured the faint-end slope of the Lya lumi- 
nosity function at redshift 5.7, using Multi-slit Narrow- 
band Spectr oscopy (MNS) to iden tify faint emission lines. 
In Paper II ( Dressier et al.|[20TT ), we presented over 200 
Lya candidates and showed that a steep faint-end slope 
(a '-^ —2) was likely. Now, we have added spectroscopic 



followup to distinguish foreground line emitters and con- 
firm LAEs. In this way, we have conclusively identi- 
fied three, and tentatively identified three further LAEs. 
With these LAEs in hand, we develop the appropriate 
framework for deriving a luminosity function from blind- 
spectroscopic data, where the positions of galaxies within 
the search slit are not always known. The resulting Lya 
LF has a faint-end slope of -l.lOtofr (-l-45to:7o) '^hen 
six (three) LAEs are used. When our spectroscopic fol- 
lowup in COSMOS and the LCIRS 15H field is complete, 
we can expect to identify around 40 LAEs, and improve 
the uncertainties on the LF measurements. 

Studies of galaxies in the post-reionization era (z ~ 
4 — 6) allow us to map the growth of young objects in an 
early universe. For the LAEs in our survey, we estimate 
SFRs of one to a few Mq yr~ ^ , and a stellar mass of at 
least 10* M0. In fact, these properties are matched to the 
masses, SFRs, and redshifts of galaxies that are thought 
to have enriched the IGM w ith metals ( Oppenheimer et 
al]|2009[ partin et al.|[2010| . 

By identifying some ot the faintest galaxies ever ob- 
served at z ~ 6, we uncover objects that are important 
for the reionization of the IGM. In |8.5[ we show that 
by extrapolating our faint-end slope measurement to low 
luminosities, it may be possible for LAEs to maintain an 
ionized IGM with our nominal assumptions about the 
escape fractions of LyC and Lya photons {fLyC — 0.1 
and fhya — 0.5) and the dumpiness of the ionized hy- 
drogen in the IGM (C = 6). Nevertheless, the allowed 
Lya luminosity density of sources brighter than 10^^ erg 
s^^ spans more than an order of magnitude. If the true 
value falls among the lower range allowed by our data, 
the escape fractions and clumping factor may need to be 
revised in order for galaxies to maintain an ionized IGM. 
Ultimately, this requirement may prove reasonable- the 
escape fractions and clumping factor are all uncertain. 
Future efforts to constrain these parameters, in parallel 
with improved measurements of galaxies will be needed if 
we are to conclusively determine, whether star-formation 
at z ~ 6 maintains the ionized IGM. 

The LAEs presented here set a benchmark for compar- 
ison to higher redshift Lya surveys that aim to measure 
the reionization of the IGM. By combinin g our line width 
measurements with those reported by Hu et al. (2010), 
we confirm that fainter LAEs have somewhat narrower 
line widths. Although the origin of this correlation is un- 
clear, it must be taken into consideration when drawing 
conclusions about a (partially) neutral IGM from Lya 
line profiles, since higher-redshift surveys will naturally 
be limited to more luminous LAEs. Furthermore, our 
measurement of the faint-end slope of the Lya LF repre- 
sents the first step towards measuring the faint-end evo- 
lution imparted by a neutral IGM. An MNS survey at 
redshift 6.5 cou ld determine if the evolution observed at 
the bright-end ( jKashikawa et al.|2006 2011 Ouchi et al. 
2010 1 persists at the tamt-end. Such observations may 



constrain models that describe the propagation of ionized 
regions in a partially neutral IGM (e.g. M iralda-Escud^ 
et al.|2000{|Ciardi fc Madau|2003||Furlanetto et al.|200' 
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In summary, we have shown that spectroscopic 
searches for LAEs are able to uncover faint galaxies that 
cannot be observed by any other means. Measuring this 
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population opens new windows to understanding the for- 
mation of young galaxies, and the enrichment and reion- 
ization of the IGM. 
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Fig. 14. — The source counts of all emission line galaxies in our survey follow a power law, and are highly complete brighter than 
F = 10~^^'^ erg cm~'^. A fit to the counts above this threshold gives a power law index of —1.02 it 0.12. Comparison of the observed 
counts to the extrapolated fit gives the survey completeness, shown on the right. 
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APPENDIX 
A. SURVEY COMPLETENESS 



To calculate the survey incompleteness for the IMACS search data, we take advantage of the fact that our emission 



line counts are well described by a power law between F = 10~^^ and 10~^^-^ erg s~^ cm~^ (see Figure 14) . This is 
easily understood, because we are summing the faint-end counts of (primarily) the constituent foreground populations- 
all of which are described by power-laws. Clearly, we can expect these counts to continue to rise to fainter emission 
line fluxes. Therefore, we predict the the number of all emission line galaxies by fitting power-law to the counts 
brighter than 10~^^'^ erg s~^ cm~^ in both the COSMOS and 15H fields. (The sensitivity is the same for both fields, 



so we include the 15H field here to improve statistics.) This fit, illustrated in Figure |14| goes as ^-i 02±o.i2^ r^j^^ 
completeness of our survey- defined as the ratio of the observed counts to those predicted by our model- is plotted 
in Figure |14| A completeness fraction slightly greater than unity naturally arises when the power-law fit falls slightly 
below the oata. 

For comparison, we show the model adopted by Paper II, which was estimated from a modification to the com- 
pleteness function presented in Paper I. The good agreement between the estimate in Paper II and the completeness 
calculated here indicates that no substantial revisions to the implications in Paper II are needed. Additionally, we 
note that the faintest bin that we consider in this paper contains many galaxies (42), so the Poisson noise on this 
calculation makes a negligible contribution to our error budget. Finally, while one might expect that the completeness 
to line emission is different for LAE candidates (without detected continuum) and foreground galaxies (with detected 
continuum), the faintest bins are dominated by LAE candidates, so it is reasonable to apply this incompleteness 
correction to the LAEs in our survey. 

B. FOREGROUND POPULATIONS 

In Paper II, we characterized the foreground Ha, [O III], and [O II] populations that comprise our sample of emission 
lines. Subtracting the foreground LFs estimated from the COSMOS field showed a significant excess of counts that 
suggested a steep rising slope {a ~ —2) for the remaining galaxies that are likely LAEs. Additionally, an angular 
correlation of the Lya candidates and [O II] emitters showed that fewer than 30% of the LAE candidates should be 
[O II]. Since even fewer should be Ha and [O III], the implication is that approximately 100 LAEs may exist in the 
combined COSMOS and 15h survey fields. In this section, wc show that the foreground number counts in COSMOS 
are indeed higher than what we adopted in Paper II, consistent with our preference for a slightly shallower faint-end 
slope for the Lya LF. 

In f|5] we identified foreground galaxies through two means: objects which are obvious foreground emitters in the 
search data, and objects which were initially LAE candidates in the search data, but are now spectroscopically or 
photometrically identified as foreground galaxies. The counts of the total foreground populations are listed in Table [3j 
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TABLE 6 

Foreground Number Counts 



Flux 


Identified Forej 


;rounds Unidentified A^pred^ 


Inferred (Low) 


Inferred (High) 


A^pred (+ lo-)'' 


-17.47 


19 


22 


23.0 


38.9 


42.0 


33.6 


-17.27 


39 


13 


35.8 


47.8 


49.9 


49.0 


-17.07 


33 


9 


30.0 


38.2 


40.8 


38.3 


-16.87 


24 


6 


22.5 


~ 30 


~ 30 


27.0 


-16.67 


14 


2 


16.5 


~ 16 


~ 16 


18.5 


-16.47 


14 





10.8 


14 


14 


11.3 


-16.27 


13 





6.2 


13 


13 


6.3 


-16.07 


6 





3.0 


6 


6 


3.0 


-15.87 


2 





1.4 


2 


2 


1.4 


-15.67 


1 





0.7 


1 


1 


0.7 


Model 


predictions of tlie 


observed foreground 


galaxy counts 


in include the incompleteness given 


in Appendix [a| 



and shown again, for reference in Table [6] where we compare these counts to the numbers predicted by the foreground 
luminosity functions used by Paper II. In each bin, the foreground galaxy counts that we have identified so far are 
similar to (often exceeding) the number predicted by foreground LPs. In addition, in f|5] we showed that the results 
of our spectroscopic followup indicate that the unidentified LAE candidates (for which counts are also listed in Tables 
[3] and |6|) are predominantly foreground galaxies. Consequently, we expect that when our spectroscopic followup is 
complete, the observed foregrounds will be higher than the model favored in Paper II. We estimate these counts in 
the same way that we calculated the inferred LAE counts in f|5j The inferred foreground number counts are given for 
two cases, with the low (high) counts referring to the case of six (three) genuine LAEs in this paper. 

Including the foreground galaxies inferred among the unidentified line emitters makes the case clear: the nominal 
"best estimate" foreground counts in Paper II were underestimated, at least for the COSMOS field. However, they 
remain within range of reasonable foreground LPs which we considered. An alternative model is shown in Table |6] 
where the faint-end slopes are all steeper by la from their adopted values (see Paper II). These counts agree more 
closely with our expectation from followup spectroscopy. However, a model where all three foreground populations 
have steeper slopes by la is contrived. It is more probable that the number counts of one of the foreground populations 
has been significantly underestimated. 

The breakdown of the spectroscopical ly identified fo reground galaxies in Table H] suggests that H/3 emitters are the 
cause of the larger foreground counts. In Dressier et al. (j2011j, this foreground population was not explicitly considered. 
It was, however, implicitly included with the [O IIIJ population. When foreground emission line galaxies are detected 
throi igh narrowband imaging, broad-band colors are typically used to identify the approximate redshift and probable 
line ( [Takahashi et al.||2007| |Ly et al.|[2007[ [Shioya et al.|[2008l ). Not only does this mean that the [O III] catalog wih 
contain some H/3 emitters; an Ha catalog may be contaminated by [S II] 6716, 6730 A emitting galaxies. In practice, 
spectroscopic followup of bright emission lines from narrowband imaging samples has shown that contamination of 
this form is rare (Ly et al. 2007 find 16 [O 111] emitters for only two H/3 emitters). Even if this contamination were 
significant, the overall [O III] -|- H/3 luminosity function should be sufficient to allow for an accounting of the total 
foreground populations at the bright end, as we did in Paper II. 

Our spectroscopic followup observations, which included known foreground objects as slit-mask filler, provide strong 
evidence that H/3 emitters in COSMOS play a more predominant role than we assumed in Pap er II . In total, we have 
identified 15 [O III] and 21 H/3 emitters. The distributions of their fluxes are shown in Figure 15 
the numbers of [O 111] emitters, an excess of H/3 emitters fainter than about 8 x 10""'^* erg s 



where, relative to 
^ is observed. 



10~"" erg s~" cm 

One way to interpret the excess of faint H/3 (relative to [O III] ) emitters is as an indication that the faint-end of the 
H/3 LP is significantly steeper than the [O III] LP. This would be the case if fainter [O III] emitters have progressively 
lower [O lll]/H/3 ratios. Since [O llI]/H/3 is strongly correlated with the R23 ({[O 111]+ [O ll]}/H;g; |Pagcl ct al. 19791 ) 
metallicity diagnostic, low [O 111] /H/3 ratios could be an indication of low oxygen abundance (when galaxies are on 
the lower-branch of R23). In prac tice, however, even the lowest metalli city galaxies observed at z < 1 have high 
[O III]/H/3 ratios (larger than three; [Salzer et aL|2005[ [Kakazu et~al]|2007| . The [O III]/H/3 ratios that we measure m 
our DEIMOS observations are similarly high, 'lliis nieans that it is ditticult to interpret the excess H/3 counts as an 
indication of low-metallicity galaxies. 

A more probable explanation for the H/3 emitters is that we hav e observed an ove rdensity at z • 
we use the zCOSMOS (bright) catalog of spectroscopic redshifts (jLilly et al. 



0.68. 



In Pigure [T5| 

|2007[ ) to show the spatial distributions 
of galaxies in the redshift windows that corresponds to H/3 and [0 lllj AA4959, 5007. Even though the volume at the 
"[O 111] redshift" is 50% larger than at the "H/3 redshift", the full zCOSMOS catalog contains 2.5 times as many 
galaxies in the latter redshift window. Within the footprint of our IMACS field (see Pigure [Tsj right), this excess 
density is amplified, with a factor of f our more galaxies in the redshift window at z ^ 0.68 than at z ^ 0.63. The 



excess of faint H/3 emitters in Pigure 15 (left) is easily explained by this overdensity, which is most apparent at 
luminosities fainter than the knee of the~H/3 LP. Therefore, the COSMOS foreground counts at faint fiuxes are indeed 
higher than the model adopted in Paper II. 

This increase over our previous estimates of the foreground populations in COSMOS does not erase our signal from 
LAEs. Although we find a slightly shallower slope (—1.7 as opposed to —2.0), these results are consistent, as a = —2.0 
falls within the 68% confidence interval that we present in this paper. Purther followup spectroscopy (including our 
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Fig. 15. — Left- Spcctroscopically confirmed H/3 emitters outnumber [O III] emitters at faint fluxes. Right- The zCOSMOS (bright) 
spectroscopic catalog (Lilly et al. 2007) contains a significantly higher density of galaxies at 2: ~ 0.68 than at 2 ~ 0.63. Objects are shown 
for the redshift windows that correspond to H/9 and [O III]AA4959, 5007 in our survey, although the latter volume is 50% larger. The 
footprint of our IMACS field is shown by the black circle. Within this region, there are four times as many galaxies at 2 ~ 0.68 than at 
2 ~ 0.63, implying a volume density six times higher. We note that the coverage of zCOSMOS is not uniform over the regions shown here, 
but the relative numbers of the green and blue points should be unaffected. 

15H field) will constrain the faint-end of the Lya LF more tightly and simultaneously measure the faint foreground 
populations. 



